Nuclear Charge Radii of 911 Li: the Influence of Halo Neutrons 

R. Sanchez, 1 W. Nortershauser, 1 * 2 G. Ewald, 1 D. Albers, 3 J. Behr, 3 P. Bricault, 3 B.A. Bushaw, 4 A. Dax/'B 

J. Dilling, 3 M. Dombsky, 3 G.W.F. Drake, 5 S. Gotte, 1 R. Kirchner, 1 H.-J. Kluge, 1 Th. Kiihl, 1 J. Lassen, 3 
CD. P. Levy, 3 M.R. Pearson, 3 E.J. Prime, 3 V. Ryjkov, 3 A. Wojtaszek/'Q Z.-C. Yan, 6 and C. Zimmermann 2 

1 Gesellschaft fur Schwerionenforschung, 64291 Darmstadt, Germany 
Departement of Physics, University Tubingen, D-72076 Tubingen, Germany 
3 TRIUMF, Vancouver, BC, Canada V6T 2 A3 
4 Pacific Northwest National Laboratory, P.O. Box 999, Richland, WA 99352, USA 
J Department of Physics, University of Windsor, Windsor, Ontario, Canada, N9B 3P4 
6 Department of Physics, University of New Brunswick, Fredericton, New Brunswick, Canada E3B 5A3 

(Dated: February 9, 2008) 

The nuclear charge radius of 11 Li has been determined for the first time by high precision laser 
spectroscopy. On-line measurements at TRIUMF-ISAC yielded a 7 Li — u Li isotope shift (IS) of 
25101.23(13) MHz for the Doppler-free 2s 2 S 1/2 -> 3s 2 S 1/2 transition. IS accuracy for all other 
bound Li isotopes was also improved. Differences from calculated mass-based IS yield values for 
change in charge radius along the isotope chain. The charge radius decreases monotonically from 
6 Li to 9 Li, and then increases from 2.217(35) fm to 2.467(37) fm for xl Li. This is compared to 
various models, and it is found that a combination of halo neutron correlation and intrinsic core 
excitation best reproduces the experimental results. 



PACS numbers: 32.10.Fn, 21.10.Ft, 27.20.+n 

For twenty years, halo nuclei with diffuse outer neu- 
tron distributions have been known to exist at the limits 
of stability for many of the lighter elements 0. The 
first discovered [2j and most renowned of these is 11 Li 
with two halo neutrons; however, details of the nuclear 
structure and halo - core interactions are still not well 
understood. Nuclear forces are not strong enough to bind 
a neutron to 9 Li, nor can they bind two neutrons into a 
dineutron. Yet adding two neutrons to 9 Li leads to a 
bound nucleus — n Li (T x / 2 = 8.4 ms), illustrating the 
importance of understanding the interaction that allow 
formation of the halo structure. Recent measurement |3j 
of the rms nuclear charge radius (r c ) for the two-neutron 
halo 6 He indicates that its halo is a dineutron "orbiting" 
the He core. The core is a strongly-bound a-particle 
and model calculations |4| estimate only a 4% increase 
in r c (a). In contrast, the 9 Li-like core of n Li is 'softer' 
and interaction between halo neutrons and core nucleons 
may significantly polarize the core. 

An indicator for an altered 9 Li core would be a change 
in proton distribution between 9 Li and 11 Li. This was 
investigated in collisions which removed a proton from 
a n Li projectile [f|, but within the rather large uncer- 
tainty, there was no clear evidence for a change in the 
deduced charge radius. Also, analysis and interpretation 
was not straightforward because of the dependence on an 
assumed nuclear model. A more sensitive approach to de- 
termine the change in r c is a measurement of the isotope 
shift in an atomic transition A finite nuclear charge 
distribution reduces electron binding energies, particu- 
larly for s-clcctrons that have probability of being inside 
the nucleus, and a change in the distribution between 
isotopes can be observed as shifts in electronic transition 



energies. In light elements, the mass-based isotope shift 
is much larger than the nuclear volume shift; for lithium, 
about 10,000 times larger. The dominant portion of the 
mass shift is change in reduced mass (normal mass shift), 
but electron correlations (specific mass shift) are also im- 
portant. Recent high-precision calculations account for 
these correlations, as well as relativistic and QED cor- 
rections [3. In this work we present the first measure- 
ment of the n Li isotope shift in the 2s — > 3s Doppler- 
free two-photon transition, as well as refined values for 
all other isotopes. These are compared with calculated 
mass shifts, yielding nuclear charge radii that are com- 
pared with various theoretical models and interpreted in 
terms of halo correlation and core polarization. 



The experiment must fulfil two conditions: measure 
the isotope shift to an uncertainty of one part in 10 5 
(10~ 10 of total transition frequency), and provide an 
overall efficiency sufficient to observe the resonances with 
production yields of ~ 10 4 n Li atoms/s. Moreover, the 
short half-life requires on-line study at the production fa- 
cility. We previously reported a technique [8j to perform 
such measurements and used it to determine the charge 
radii of 8 - 9 Li produced at the GSI-UNILAC. For the ex- 
periments reported here, the apparatus was moved to the 
TRIUMF-ISAC facility in Vancouver, Canada where n Li 
is produced by a 40 pA, 500 MeV proton beam imping- 
ing on a tantalum target. n Li + ions extracted from the 
target (typ. 30,000/s) are implanted in a hot carbon foil 
where they are neutralized and released as atoms into 
the low-field source region of a quadrupole mass spec- 
trometer (QMS). The neutral atoms are re- ionized via 
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FIG. 1: Resonances in the 2s — > 3s transition of 11 Li as a 
function of the beat frequency between the titanium-sapphire 
laser and the reference diode laser. Error bars are simple 
counting statistics on the number of observed ion counts. 



TABLE I: Isotope shifts measured at TRIUMF (this work) 
and GSI 8] [avg = weighted mean] compared with theoretical 
mass shifts for 7 Li — A Li in the 2s 2 <Si/2 — * 3s 2 Si/2 transi- 
tion. Uncertainties for r c are dominated by uncertainty in the 
reference radius r c ( 7 Li) = 2.39(3) fm [Til ]. 



Isotope Isotope Shift, kHz Mass Shift, kHz 



fm 



'Li 



s Li 



'Li 



this 
GSI 
avg 
this 
GSI 
avg 
this 
GSI 
avg 
this 



-11453 984(20) 
-11453 950(130) 
-11453 983(20) 
8 635 781(46) 
8 635 790(150) 
8 635 782(44) 
15 333 279(40) 
15 333140(180) 
15 333 272(39) 
25 101 226(125) a 



-11453 010(56) 2.517(30) 



8 635 113(42) 2.299(32) 



15 332 025(75) 
25 101812(123) 



2.217(35) 
2.467(37) 



kHz statistical + 57 kHz systematic from AC-Stark shift 



doubly-resonant four-photon ionization: 

r, 2x735nm „ decay 610nm „ , 610,735nm i , » 

2s > 3s > 2p > 3d ■ > LL , (1) 

30ns 

with excitation taking place at the center of a doubly- 
resonant optical enhancement cavity (^100 x) built 
around the QMS source region. The titanium-sapphire 
laser that excites the 2s — > 3s two-photon transition is 
beat-frequency servo-locked to an I2 hyperfine line stabi- 
lized diode laser. As previously described 0, measure- 
ments on 11 Li (and the other isotopes) were interspersed 
with measurements on 6 Li which served as the experi- 
mental reference, and measured optical powers were used 
to correct for calibrated AC-Stark shifts. 

Figure^ shows a typical n Li spectrum. Twenty-four 
such spectra were obtained over six days of beam time. 
With nuclear spin / = 3/2, the 2 S1/2 states have F = 1, 2 
hyperfine components, which obey the two-photon se- 
lection rule AF = for an s — ► s transition. All Li 
isotopes have nuclear spin and exhibit similar doublets: 
Isotope shifts are taken with respect to center-of-gravities 
of the two hyperfine lines for each isotope. Results for 
all isotopes, relative to 7 Li, are given in Tabled Val- 
ues for 6:8,9 Li are in good agreement with our previous 
measurements 0, but with improved precision. The 6 Li 
isotope shift was also determined earlier with a different 
technique as -11453 734(30) kHz |9J; this is significantly 
different from our current measurements (~ 5 times 
the combined uncertainties), and is attributed to un- 
accounted systematic errors in the prior interferometric 
measurements 9] , as compared to the current frequency- 
based measurements. The IS for the halo nucleus n Li is 
a first-time measurement. 

Successful determination of changes in r c from the iso- 
tope shift measurements depends critically on the com- 
bined accuracy of theory and experiment. On the the- 



oretical side, the quantum mechanical many-body prob- 
lem must be solved to high accuracy in the nonrelativis- 
tic limit, and then the effects of relativity and quantum 
electrodynamics are included with perturbation theory. 
In the past, theoretical results with laser-spectroscopic 
accuracy were not available for atoms more complicated 
than helium, even in the nonrelativistic limit. This prob- 
lem is now solved by variational methods involving cor- 
related basis sets with multiple distance scales 6J. The 
resulting electron wave functions are used to calculate the 
various contributions to the mass shift, listed for 7,11 Li 
in TablelTTI A recent first calculation of the mass 
polarization correction to the Bethe logarithm part of 
the electron self-energy has significantly reduced uncer- 
tainty in the QED contribution; overall calculation un- 
certainty is now limited by the relativistic recoil term of 
order a 2 (/i/M). 

The total in TablelTTI is the calculated mass-based com- 
ponent of the isotope shift; corresponding shifts for all 
isotopes are obtained directly from coefficients given in 
Table III of Ref. [7j and are listed in Tabled] Differences 
from measured isotope shifts are then attributed to the 
nuclear volume effect and are related to r c of the two 



TABLE II: Contributions to the 7 Li — 11 Li mass shift in 
the 2s 2 Si/2 — » 3s 2 Si/ 2 transition, excluding nuclear size 
effects. /i/M is the ratio of the reduced mass and the atomic 
mass. Uncertainty in the nonrelativistic (fi/M) term is from 
uncertainty in the 11 Li mass [To| . while limiting uncertainty 
in the relativistic and QED terms is computational. 



Contribution (order) 


kHz 


Nonrelativistic (fj,/M) 
Nonrelativistic (/x/M) 2 
Relativistic a 2 (fi/M) 


25 104 483 (20) 


-2 968 (0) 


417(121) 


QED a 3 ((jL/M) 


-120 (6) 


Total 


25 101812 (123) 



3 



2.7 
2.6 
2.5 
| 2.4 
2.3 
2.2 
2.1 



— 1 


1 


1 1 


1 


1 — 

o 










V 














g... 


A \S/ 




T 












1 

6 


i 

7 


i i 

8 9 


1 

10 


1 

11 






Li Isotope 







FIG. 2: (Color online) Experimental charge radii of 
lithium isotopes (•) compared with theoretical predictions: 
A: Greens- Function Monte Carlo Calculations |jl Hsl |. V: 
Stochastic Variational Multi-Cluster Model EE E3 (▼: as- 
suming a frozen 9 Li core), ffi: Fermionic Molecular Dynamics 
[Tsll . o: Dynamic Correlation Model fTgj . □ and o: ab-initio 
No-Core Shell Model HHH. 
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where Ze is the nuclear charge and (5(ri)) are expec- 
tation values for electron density at the nucleus in the 
respective states 0. 

Optical isotope shift measurements provide only the 
change in the rms nuclear charge radius between two iso- 
topes. Absolute charge radii r c must be referenced to at 
least one isotope that is determined with a different tech- 
nique. For the stable 6 ' 7 Li isotopes, rv have been deter- 
mined by elastic electron scattering |ll| . from which we 
use r c ( 7 Li) = 2.39(3) fm as a reference radius. This and 
the measured 6,7 Li isotope shift yields r c ( 6 Li) = 2.52(3) 
fm, in good agreement with the electron scattering result 
of 2.55(4) fm 0. Combining measured isotope shifts, 
calculated mass shifts, and the 7 Li reference radius yields 
r c for the other isotopes, as given in the last column of 
Tabled 

The derived nuclear charge radii are shown as filled 
circles in Fig. [21 while r c decreases continuously from 
6 Li to 9 Li, there is a large increase from 9 Li to n Li. 
The significance of these results becomes evident when 
compared with predictions from different nuclear models, 
also shown in Fig.[2J Models using point-proton radii 
r pp are converted to nuclear charge radii r c by folding in 
proton \\2\ and neutron [l3T | mean-square charge radii: 



iV 



(rl) = (r 2 pP ) + « + V W + 



3h 2 



where the last term is the Darwin-Foldy correction for 
"Zitterbewegung" of the proton [Tij ]. 

Neither conventional shell model nor self-consistent 
Hartree-Fock calculations have correctly reproduced 
halo-specific anomalous properties of light nuclei close to 
the neutron drip-line. Early models for n Li only treated 
its three-body character, without considering possible 
polarization of the 9 Li core j2^; thus, change in nuclear 
charge radius could only be caused by correlation of the 
two halo neutrons. If they spend most of their time on 
the same side of the core, the center-of-mass (CM) is 
clearly different from the core center, the 9 Li core orbits 
the CM, and the averaged charge distribution is diffused. 
Forssen et al. |2^| constructed corresponding wave func- 
tions for 11 Li to obtain an analytical model for electro- 
magnetic dissociation of halo nuclei. The n Li rms matter 
radius of 3.55 fm was adjusted to be in good agreement 
with experiment j2^|; the predicted CM - core distance 
Rcm ranged from 0.8 fm to 1.08 fm. The approximation 



Li) 



~ 2 ( 9 Li)] 



1/2 



2.40(6) fm is in 



4m^c 2 



(3) 



reasonable agreement with, but slightly lower than, our 
experimental result. However, information available on 
the binary neutron- 9 Li (core) interaction is insufficient 
for these calculations to yield structural details on n Li, 
nor do they make predictions for changes in r c between 
the non-halo nuclei. 

The dynamic correlation model (DCM) is a more ad- 
vanced scheme that starts from shell model states, and 
then introduces neutron-core interaction with a two-body 
potential 0,H(|. This leads to an admixture of virtually 
excited single-particle states from the core. For 11 Li, ex- 
cited bound and continuum states of 9 Li up to 50 MeV 
were included in the analysis. Charge radii calculated 
(o in Fig.|2l for 6 ' 7 Li agree well with our measurements, 
and while those for 9 ' n Li are clearly overestimated, the 
increase from 9 Li to n Li is correctly reflected. 

More sophisticated nuclear models treat interactions 
between individual nucleons using realistic nucleon- 
nucleon (NN) and three-nucleon (NNN) interactions. NN 
potentials are usually based on the multi-energy partial- 
wave analysis of elastic NN scattering data produced by 
the Nijmegen group [23] in 1993, while the NNN interac- 
tion parameters are adjusted to fit the binding energies of 
light nuclei. Greens-Function Monte-Carlo (GFMC) cal- 
culations 0, Il5l ] , the most fundamental description avail- 
able for light nuclei, have been completed for most nuclei 
with mass numbers A < 12. Results for the isotopes 
6 ' 7 ' 8,9 Li are shown (A) in Fig. [21 and are in good agree- 
ment with the experimental results. The general trend is 
reproduced, but thus far the model has not been able to 
reproduce the 11 Li binding energy. 

No-Core Shell Model (NCSM) calculations have been 
performed using realistic NN potentials. Earlier calcula- 
tions [2l| (o) for 7 > 8 > 9 ' 1:L Li treated three-body interactions 
as an effective phcnomcnological potential, while recent 



4 



work [2(| (□) included microscopic three-body potentials 
and was applied to 6 ' 7 Li. As seen in Fig.0 neither the 
absolute charge radii, nor the trend along the isotopic 
chain are in agreement with our results. 

The Fermionic Molecular Dynamics (FMD) model 
uses Gaussian wave packets for individual nucleons. The 
Winteraction is derived from the Argonne V18 inter- 
action, treating short-range correlations explicitly with 
a unitary operator. Predictions of the model [18j (©) 
are in good agreement with experiment for 6,7,8 ' 9 Li, but, 
like GFMC, the halo structure of n Li has not yet been 
successfully modelled. 

Calculations that consider interactions between all 
individual nucleons quickly become very complex and 
time consuming with increasing nucleon number. Clus- 
ter models like the stochastic variational multi-cluster 
(SVMC) calculations of Varga et al. 0,0] freeze some 
parts of the model space and allow focus on those degrees 
of freedom thought to be most relevant to the physical be- 
haviour of a given nucleus. To a large extent, this cluster 
structure can also be identified in FMD calculations. The 
building blocks in the SVMC model are the nucleons p 
and n, the a-particle, and the tritium nucleus t. The nu- 
clei a and t are not treated as structureless particles; their 
wave functions are constructed on the nucleonic level and 
only nucleon motion within the clusters is approximated 
by simple shell-model configurations. The many-body 
state then describes the correlated relative motion of the 
different clusters in a fully anti-symmctrized wave func- 
tion that obeys the Pauli principle and thus also accounts 
for correlated motion of the halo neutrons. The nucleon- 
nucleon interactions are chosen to reproduce, e.g., phase 
shifts in NN, aN and aa scattering, and deuteron size 
and binding energy. Additional effective nucleon-nucleon 
interactions are included to account for three-nucleon in- 
teractions. This model clearly shows the best agreement 
with our experiment (V in Fig.EJ). Calculations for n Li 
were performed both with and without possible excita- 
tions of the 9 Li core by the halo neutrons. Including 
these intrinsic excitations results in z'c ( 11 Li) = 2.52 fm, in 
good agreement with experiment, while neglecting them 
results in the much smaller value r c ( 11 Li) = 2.28 fm (T in 
Fig.|2J. Thus, within the framework of SVMC, neutron 
correlations alone cannot reproduce the large change in 
r c between 9 Li and n Li observed in the experiment. The 
calculations rather indicate that the core is indeed per- 
turbed and that this perturbation accounts for most of 
the charge radius increase. It will be interesting to see 
whether the model is also able to describe correlations in 
the momentum distributions of breakup fragments |29| . 
We also note that while the SVMC model clearly shows 
the best agreement with our measured nuclear charge 
radii, it's predictions for nuclear electromagnetic mo- 
ments ^3 s tiU have significant discrepancies from ex- 
perimental values, indicating that further work is still 
needed. 
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